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The osmotic coefficient of B-DNA in water may, in moderately dilute solutions, deviate as much 
as 100 % from predictions based on a simple 'counterion condensation' theory. We determine the 
results for osmotic properties via a cell model description of the ionic atmosphere near a cylindrical 
polyelectrolyte. The cell model predictions for the osmotic properties disagree with predictions based 
on simple condensation theory, but are in surprisingly good harmony with experimental findings. We 
argue that the neglect of finite-radius effects makes simple condensation theory inapplicable at all 
but impractically low polyelectrolyte concentrations and, unable to reproduce osmotic properties of 
polyelectrolytes such as DNA. 
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Ionic screening of charge interactions remains one of the most vigorously discussed properties 
of polyelectrolyte solutions ^. Recently measured osmotic coefficients of B-DNA in dilute 
aqueous solutions 0, show factor-of-two deviations from predictions based on the influential 
counterion condensation theory of Oosawa and Manning 0, ^, ^]. We shall argue that an 
earlier theory pioneered by Lifson and Katchalsky based on the cell model formulation 

of the full nonlinear Poisson-Boltzmann equation , provides a more successful starting point 
than simple condensation theory for examining not only the osmotic properties of cylindrical 
polyelectrolytes in low salt conditions but also the local organization of the ionic atmosphere 
around them. By explicitly including the finite size of the cylinder, the cell model description 
fundamentally disagrees with the line charge picture of the simplest condensation theory. It is 
well known || that in the absence of added salt, the simple condensation theory is a special 
limiting case of the cell model only when the cylinder radius is sent to zero or the cell radius is 
sent to infinity. We find here that neglect of finite radius is the main reason for the failure to 
describe DNA osmotic properties. 

The osmotic pressure of a polyelectrolyte (e.g., DNA) in an aqueous solution is obtained 
from the equation of state that recognizes the degrees of freedom of both the counterion and 
the macroion. The Oosawa-Manning limit |^ which decouples the counterion atmosphere 
from the density of the macroion can thus only make sense at effectively infinite dilution. 
Osmotic pressure experiments however are performed at finite macroion concentrations. In 
these experiments ^ the osmotic coefficient of DNA in 2 and 10 mM salt solutions is seen to 
be twice that predicted by the counterion condensation theory, while the Lifson - Katchalsky cell 
model accurately predicts its magnitude. The experimental osmotic coefficient varies weakly 
with the macroion concentration in the general direction predicted by the cell model. 



The cell model considered here involves a rigid (hollow or solid), charged cylindrical polymer of 
radius a, coaxially enclosed in a cylindrical (Wigner-Seitz-like) cell of radius Rq, corresponding 
to the total system volume per polymer length, see Fig. 1. The cell acts as a neutralization 
volume for the counterions; consequently the electric field vanishes at the cell wall. Counterions 
organize within the cell according to the nonlinear Poisson-Boltzmann equation for the double 
layer electric potential, u. In units of ksT/e, in the absence of added salt, 

r dr dr 2 

Here denotes a formal 'screening length' to be determined later by the average density of 
counterions. The subscript * stands for the electrostatic potential inside and outside of the 
cylinder The potentials in turn determine the charge densities, 

n,(r) =n,,oe-"*('-) . (2) 

Here n^^ = nl/^SnlBj), and Ibj = /{eksT) is the Bjerrum length. 

Because of the major and minor grooves in DNA 0, it is desireable to treat the cylinder at 
least as partly hollow and hence to give solvent and ions access to the space within the grooves. 
We shall compute here, for simplicity, the results for cylinders that are either solid or completely 
hollow to make the maximum possible 'non-specific' accumulation of countercharge near the 
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cylinder. Counterion accumulation is determined by the solutions to the PB equation for which 
the density variation across the boundary at r = a is continuous. Uo{r > a) and Ui{r < a) are 
both expressed with respect to a zero at the cell wall (r = Rq) and /tf = = 8TrlBjn{Ro), 
where n(i?o) is the density of counterions at the cell wall. 

For a < r < Ro solution of the Poisson - Boltzmann equation yields |^ 
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n„(r)=lnl 1^ cos2(2 ln(-^)) 1 , (3) 
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[0, |] and for r < a one finds 

Uiir) =Mo + 21n(l + cr2) . (4) 

The integration constants z, Rm, k, c, and uq are obtained from boundary conditions. In 
addition to the requirement of a continuous variation of the potential at r = a, these conditions 
include dui{r)/dr\Q = 0, duo/dr\R^^ = 0, and {dUo{r)/dr — duiir) / dr)\a = 2Q/a. For DNA 
the dimensionless linear charge density Q = iBj/lpoA — 4.353 is determined by the charge 
separation Ipoi — 1.7A and the Bjerrum length Ibj — 7.14A (as in water at room temperature). 
In practice one solves the problem of fixing the parameters by an iteration that starts with a 
trial partitioning of counterions inside and outside the cylinder that is subsequently refined 
until the boundary conditions are exactly satisfied. 



The osmotic properties of charged cylinders are encoded in the osmotic pressure Tiosm = 
kBTn{Ro) on the cell wall, which, of course, also codifies the equation of state for the macroion 
at the density set by Ro- (Because the electric field vanishes at the cell wall, the only contribu- 
tion to the Maxwell stress tensor comes from the osmotic pressure.) The osmotic coefficient is 
defined as the ratio of the actual osmotic pressure to the osmotic pressure of a hypothetical gas 
of uniformly distributed counterions or, equivalently, the number density of ions n{Ro) divided 
by the total density npo4 = 1/ {Ipoat^Rq)- ■ 

Fig. 1. shows the results of a simple numerical calculation of the osmotic coefficient of 
solutions of rigid hollow cylinders as a function of the molar concentration, when the radius 
of the cylinder a = lOA, the dimensionless (nominal) line charge density Q = 4.353, and the 
molar weight M is chosen as for DNA. 

(i) The calculated osmotic coefficients vary slightly with concentration. The measured 
osmotic coefficients |^ vary even more slowly. However, simple condensation theory predicts 
no variation at all. 

(ii) In the region of significant experimental interest, i.e., for concentrations of DNA phos- 
phates from 0.1 to 0.5 M, the calculated osmotic coefficients are close to 0.3 in reasonably 
close agreement with the experimental value 0.24 Simple condensation theory predicts 
= 1/(2Q) ~0.11. 

(in) As expected in the limit a/i?o — ^ 0, our formulated osmotic coefficients converge to 
the line - charge result 0(a/i?o) — ^ 1/(2Q); it is the significant deviations that occur away 
from this infinite dilution limit that concern us. 
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Specifically, consider the asymptotic expression for the osmotic coefficient in the limit where 
a/i?o is small (the line-charge limit, or the infinite-dilution limit): The outer solution with 



give us z, as follows: 



Qo = I - ztan ( z In(^) 
V 

(R 
z ln(— — ) 
-'I'm 

k'^RI = 4(1 + z'^ 



, a , arctan '^° ) — arctan i-^ 

In — = "^—^ \iL 

^Ro^ z 



(5) 

(6) 
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When a/Ro is small, z — > 0, the arctan(-)'s may then be replaced by, respectively, —n/2 and 
n/2 so that z has a weak logarithmic dependence on a/Ro: 
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The corresponding asymptotic expression for the osmotic coefficient is 
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where G{a/ Rq) is a geometric factor, which for hollow cylinders assumes the value 1. Therefore, 
for small values of a/Ro, Eq. (H) gives 
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In other words, the osmotic coefficient develops a weak logarithmic concentration dependence. 
(If the cylinders were not completely hollow an extra concentration dependence would appear at 
high concentrations. In this case the geometric factor G{a/ Rq) would deviate from unity.) Any 
concentration dependence of the osmotic coefficient is incompatible with simple condensation 
theory. Only as a/i?o — ^ does the osmotic coefficient approach a concentration- independent 
limit l/{2Qo). 



Fig. 2 shows the osmotic pressure iTosm versus the molar concentration cdna of DNA or 
equivalently versus the radius of the cell Rq. Experimental data are available at low, 0.1 — 0.5 M 

and high concentratios, 1 — 2 M []lT|. In both ranges we find that calculated osmotic pressures 
are remarkably close to capturing the magnitude and variation of the measured pressures. In 
the low-concentration range and 2 mM salt the calculated pressures are slightly too large which 
might in part reflect a weak contribution from finite salt concentration that can be dealt with 
on the basis of a simple Donnan equilibrium picture P, O, 13]. This approach however fails 



completely at higher salts, e.g. 10 mM, where the osmotic pressure is substabtially lower than 
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the values calulated from the Donnan equihbrium. Though one could make the calculated 
values of osmotic pressure in the intermediate regime of DNA densities to be even closer to the 
data, by choosing a somewhat smaller value for the DNA radius, it turns out that in this case 
one would loose the relatively good agreement in the regime of very large DNA concentrations. 
There appears to be no simple way of adjusting the DNA parameters to get a good quantitative 



fit of the calculated osmotic pressure with experiments. As pointed out many times before |I4 
it appears yet again that non-electrostatic interactions at very large DNA concentrations make 
a significant contribution to the overall osmotic pressure in the system. 

In the high-concentration range, the calculated variation of the pressure with concentration 
is clearly slower than what one observes for the experimental data. The (small) difference 



may be of non-electrostatic origin |I^ or reflect charge-discreteness effects discussed in Ref. 



13]: salt effects are likely to be unimportant in this concentration range. The renormalization 



effects due to chain conformational fluctuations, discussed in Ref.lTT|, lead to predictions of 



decreasing rates of change of the pressure with concentration. Recent work on stretching of 
DNA at various ionic conditions |jl5| also suggests that there might be an additional strong 



coupling between DNA elasticity and electrostatics. Though the details of this coupling are 



only beginning to be elucidated [16] it is concievable that local deformations of DNA would 
change the countercharge distributions at low salt conditions and thus effect also the osmotic 
coefficient. 



To what extent do the above observations elucidate 'counterion condensation, 'simple' 0, Q or 
'extended' [|, |l8Hl|] ? 



In the 'simple condensation' picture |Q, § highly charged and rigid cylindrical macromolecules 
are portrayed as line charges with explicit neglect of finite macromolecular radius, the at- 
mosphere of counterions is divided into a condensed and osmotically inactive fraction which 
redefines the line charge density, and an 'unbound', osmotically active fraction. This 'two- 
phase coexistence' is established whenever the line-charge density (or, equivalently, the elec- 
tric field at the surface of the macromolecule) becomes large enough, quantitatively whenever 
Q = ^Bj/^POA > 1; For Q > 1, condensation will bring down the effective (dimensionless) line- 
charge density from Q to 1; a fraction fc = {Q — l)/Q = 1 — 1/Q of the counterions condense. 
The remaining fraction f = 1/Q of the charge remains unbound. If the unbound fraction is 
modeled as a polarized Debye-Huckel gas of counterions, this gas contributes to the osmotic 
coefficient (p, as (pM = / (kBTng) = f/2 = 1/(2(5). For E)NA 0m — 0.11 while a significant 
fraction / ~ 0.75 of the countercharge is 'bound'. 

In the line-charge limit ^ there is no disagreement at all between cell model predic- 
tions and predictions based on 'simple' condensation theory. In the cell model ^ an effective 
'condensation range' defines a Manning radius Rm such that the fraction of the countercharge 
contained within a shell a < r < Rm is / = 1 — l/Qo, precisely the fraction of charges that are 
predicted to condense according to the Oosawa - Manning model. Using Eqs. (||)-(|^), in the 
small a/Ro limit, Rm — -Roe~2i ~ y/R^a. That is in this limit the Manning radius disappears; 
the effects of the fraction / = 1 — l/Qo — ^ 1 — 1/Q of ionic atmosphere can be absorbed 
in a redefinition of the line charge, as in 'simple' condensation theory. However when a > 
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the Manning radius Rm is finite and even diverges as one approaches the infinite-dilution hmit 
-Ro — ^ oo. 

Is it possible to reformulate or 'extend' condensation theory to account for finite macro- 
molecular radius while retaining the essential idea of two-phase coexistence between free and 
bound fractions? If not, the organization of the counterion atmosphere around cylindrical poly- 
electrolytes and the electrostatic contribution to solution properties of polyelectrolytes, must be 
understood on the basis of a Poisson-Boltzmann or a more advanced double layer description. 
To the extent that experiments like deal with a finite concentration of macroions that con- 
tribute in and of themselves to the equation of state, i.e., to the osmotic pressure and osmotic 
coefficient, the simple condenstaion picture appears to be too drastic an idealization. 



There have been many elegant attempts to extend condensation theory [Rl 17, nR, 19 . In 



essence these are all versions of the Poisson Boltzmann equation solved under various conditions 
including finite radius, explicit ion simulation, condensed layers, etc. These models create a 
distinction between bound and unbound ions or speak of condensation shells within the context 
of the continuous distributions predicted. 

(i) Many have been tempted to identify, via the cell model analysis, the fraction / = 
1 — 1/Qo of counterion residing in the shell a < r < Rm with the condensed fraction, and 
so implicitly to accept that the condensed layer is a spatially extended object with a peculiar 
sensitivity to changes in, e.g., concentration of macroions. Generalizations of this approach 
have been considered where one identifies a length scale similar to Rm and the accompanying 
'condensation shell' from an analysis of the counterion distribution function versus linear or 



logaritmic |T8|, IT9[ radial distance. They provide operational definitions of a 'condensed layer' 



but these definitions are almost a matter of terminology. 

Does such 'condensation' add anything essentially new either to P-B theory or to simple 
line-charge condensation? 

None of these approaches makes it completely clear that there is a 'Debye-Hiickel' cloud of 
counterions in the outer shell Rm < r < Rq, although at some distance from the polyelectrolyte 
potentials must weaken to a degree that a D-H description becomes accurate . 

(ii) It has been argued recently ||19|], that the fact that in the infinite-dilution limit the 
contact density n(a) reaches a constant limiting value {l/ebna^){Q — 1)^/(2(5) leads one to 
expect the 'existence of a close layer that cannot be diluted away', in other words a condensed 
layer. This would point to the conclusion that counterion condensation has a reality going 
beyond that of being a mere heuristic tool. This is perhaps true but does not really prove 
the point, in the sense that one can not even in this case distinguish 'condensation' from mere 
double layer properties. 

Comparison of the measured osmotic pressures for DNA at various ionic conditions suggests 
that at low salt concentrations the 'counterion condensation' picture does not capture the 
main features of the data. The discrepancy experimental and theoretical values of the osmotic 
coefficient, amounting to a factor of two can not be ignored. The cell-model in which the 
countercharge distribution is governed by the Poisson - Boltzmann equation and where the 
macroion is modelled as a cylinder of a finite radius produces a much better fit to data. This 
almost quantitative correspondence between the data and the calculations brings back into 
focus the counterion atmosphere without any need to invoke a 'condensation' behavior. 
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Figure 1: Osmotic coefficient, = n{RQ)/npoi vs. mean counterion concentration npoi 
{cdna in molar units) and the geometry of the cell model. n(i?o) is the counterion density at 
the cell wall. The charge density on the cylinder corresponds to one elementary charge per 1.7A 
contour length (or Q=4.353). Osmotic coefficients have been calculated for a = lOA. We have 
considered the consequences of allowing the counterions to enter the interior of the cylinder 
(dashed line), or of preventing the counterions from entering the cylinder at all (solid line). 
The osmotic pressure in the limit of vanishing a (Manning condensation limit) is shown by the 
dashed line. The calculated osmotic coefficients increase as cdna increases. The experimental 
results for DNA 0] (■) in the salt independent regime of DNA concentrations show the same 
trend, but the variation is slower than for any of the calculated data. The osmotic coefficient 
is somewhat sensitive to whether the cylinder is hollow (dotted lines) or solid (solid lines). 
Expectably when the cylinder is hollow, the osmotic coefficient is smaller. The hollow/solid 
difference increases when Rq — > a. In the dilute, a/Ro — > 0, limit the computed osmotic 
coefficient reaches the constant limiting value of 1/{2Q) ~ .11 of counterion condensation 
theory. 

Figure 2: Osmotic pressure Tiosm versus the molar concentration of DN,A cdna (upper graph) 
or the radius of the cell Rq (lower graph). The charge density on the cylinder with radius 
a = lOA corresponds to one elementary charge per 1.7 A contour length (or Q = 4.353). Hollow 
cylinder, dotted line; solid cylinder, solid line. The predicted osmotic pressure in the limit of 
vanishing a (counterion condensation limit) is shown with a dashed line. Experimental data 
are represented with symbols: ▼ osmotic pressure data at 2 mM salt 0], ■ osmotic pressure 
data at DNA concentrations where salt does not matter, • osmotic pressure at very large DNA 
concentrations where again there are no salt concentration effects |jll|. The experimentally 
determined pressures for very small values of cdna can not be explained by the electrostatic 
cell model studied here. For cdna in the range 0.1 — 0.5M, the calculated pressures deter- 
mined from the electrostatic cell model seems to capture the trend in the experimental data 
but the calculated pressures are slightly larger than the experimentally determined pressures, 
irrespective of solid-cylinder (solid line) or hollow-cylinder (dashed line) assumptions. When 
Cdna approaches 1 M , we find that the calculated pressures and the experimentally deter- 
mined pressures are not very different in magnitude, but the calculated variation of pressure 
with concentration does not fit the trend in the experimental data. 
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